KA-TP-28-2012, MAN/HEP/2012/09, MZ-TH/12-27, 
CERN-PH-TH/2012-194, IPPP/12/50, DCPT/12/100, 
FTUV-12-0709, LPN12-078, TTK-12-32 



Release Note - Vbfnlo-2.6.0 

K. Arnold^, J. Bellm^, G. Bozzi^, F. Campanario^, C. Englert^, B. Feigl^, J. Frank^, 
T. Figy^ B. Jager^ M. Kerner\ M. Kubocz^ C. Oleari^ S. Palmer\ M. Rauch^ 
H. Rzehak*, F. Schissler^, O. Schlimpert^, M. Spannowsky^, D. Zeppenfeld^ 

^ Institut fiir Theoretische Physik, Universitat Karlsruhe, Karlsruhe Institute of 

Technology, 76128 Karlsruhe, Germany 
^ Dipartimento di Fisica, Universita di Milano and INFN, 20133 Milano, Italy 
^ Institute for Particle Physics Phenomenology, University of Durham, Durham, 

DHl 3LE, United Kingdom 
^ School of Physics and Astronomy, The University of Manchester, Manchester, 

M13 9PL, United Kingdom 
^ Institut fiir Physik (THEP), Johannes-Gutenberg- Universitat, 55099 Mainz, Germany 
^ Institut fiir Theoretische Teilchenphysik und Kosmologie, RWTH Aachen University, 

52056 Aachen, Germany 
^ Dipartimento di Fisica, Universita di Milano-Bicocca and INFN, Sezione di 
Milano-Bicocca, 20126 Milano, Italy 
^ CERN, CH-1211 Geneva 23, Switzerland (on leave from: Physikalisches Institut 
Albert-Ludwigs-Universitat Freiburg, Hermann-Herder-Str. 3, 79104 Freiburg im 

Breisgau, Germany) 

Abstract 

Vbfnlo is a flexible parton level Monte Carlo program for the simulation of 
vector boson fusion (VBF), double and triple vector boson (plus jet) production in 
hadronic collisions at next-to-leading order (NLO) in the strong coupling constant, 
as well as Higgs boson plus two jet production via gluon fusion at the one-loop 
level. This note briefly describes the main additional features and processes that 
have been added in the new release - Vbfnlo VERSION 2.6.0. At NLO QCD 
diboson production {W'j, WZ, ZZ, Z7 and 77), same-sign W pair production via 
vector boson fusion and the process Wjjj have been implemented (for which one- 
loop tensor integrals up to six-point functions are included). In addition, gluon 
induced diboson production can be studied separately at the leading order (one- 
loop) level. The diboson processes WW, WZ and Wj can be run with anomalous 
gauge boson couplings, and anomalous couplings between a Higgs and a pair of 
gauge bosons is included in WW, ZZ, Zj and 77 diboson production. The code 
has also been extended to include anomalous gauge boson couplings for single vector 
boson production via VBF, and a spin-2 model has been implemented for diboson 
pair production via vector boson fusion. 
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1 Introduction 



Vbfnlo |T|[2] is a flexible Monte Carlo (MC) program for vector boson fusion (VBF), 
double and triple vector boson (plus jet) production processes at NLO QCD accuracy. 
The electroweak corrections to Higgs boson production via VBF have been included. In 
addition, the simulation of CP-even and CP-odd Higgs boson production in gluon fusion, 
associated with two additional jets, is implemented at the (one-loop) LO QCD level. 
Vbfnlo can be run in the MSSM, and anomalous couplings of the Higgs boson and 
gauge bosons have been implemented for certain processes. Additionally, two Higgsless 
extra dimension models are included - the Warped Higgsless scenario and a Three-Site 
Higgsless Model - for selected processes. 

Further information, and the latest version of the code, can be found on the Vbfnlo 
webpage 

'http: //www- itp. particle .uni-karlsruhe .de/vbfnlo/ 

A complete process list is given in Appendix A. 

2 New processes 

The latest version of Vbfnlo has several new processes implemented at NLO QCD. 
2.1 Same sign W pair production via VBF 

Same sign W pair production with two jets via vector boson fusion js] has been included in 
Vbfnlo 2.6.0. This process is potentially sensitive to new physics signals and, as it gives 
rise to same sign dilepton final states, it is also a background to new physics scenarios. 
To distinguish potential signatures of physics beyond the Standard Model from the effect 
of higher order corrections, precise theoretical predictions are essential. Although the 
QCD corrections to the integrated cross sections were found to be relatively small, their 
effect on several distributions is appreciable. When NLO QCD corrections are taken into 
account, the residual scale uncertainties are at the 2.5% level. The new process IDs are 
given in Table [Tj 



ProcId 


Process 


250 


pp W+W^ 


' 33 - 


-> tl^i>^l2^i>^33 


260 


pp — ^ W^W' 


' 33 - 





Table 1: New process IDs for diboson + 2 jet production via vector boson fusion at NLO 
QCD accuracy. 
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2.2 Diboson production and gluon-induced contributions 



A good understanding of diboson processes at the LHC is essential, as not only do they 
allow the study of the Standard Model's gauge structure, but they also provide a back- 
ground to Higgs boson and new physics searches. Anomalous triple gauge boson couplings 
are included in the WZ and W'j processed (For contributions from anomalous HVV 
couplings in gluon induced processes - Fig. flp see below.) The new processes are included 
at NLO QCD accuracy under the process IDs shown in Table |2} 



ProcId 


Process 


BSM 


310 


(-) 
PP - 

(-) 




Z 




|> anomalous VVV couphngs 


320 


PP - 


w- 


Z 




330 


(-) 
PP - 


-> zz 






anomalous HVV couplings 


340 


(-) 
PP - 

(-) 




7- 




!> anomalous VVV couplings 


350 


PP - 


w- 


7 - 




360 


{-) 
PP - 

(-) 


Z'j 


— )- 




j> anomalous HVV couplings 


370 


PP - 


^ 77 







Table 2: Process IDs for the new diboson production processes at NLO QCD accuracy. 



In addition to the NLO QCD corrections, the gluon-induced fermionic loop processes 
can be included in those processes with neutral final states (including WW production, 
which was already implemented in Vbfnlo, with process ID 300). Both continuum 
production via box diagrams as well as production via an s-channel Higgs boson resonance 
(shown in Fig. [T]) are available, with full interference effects. 





Figure 1: Gluon-induced contributions to diboson production. 

The vbfnlo.dat flag FERMIONLOOP controls these contributions: 

• switches off these processes 

• 1 includes only the box contribution 

• 2 includes only diagrams via an s-channel Higgs resonance 

• 3 includes both contributions including interference effects. 



^As in the rest of Vbfnlo, no neutral triple gauge couplings are included. 



Although these contributions are formally of higher order, their effect can still be signif- 
icant. Anomalous couplings between a Higgs boson and a pair of gauge bosons can be 
included in those processes with neutral final states. The gluon-induced loop diagrams 
can be studied separately in the new release of Vbfnlo using the process IDs of Table |3] 
and the executable ggf lo. 



ProcId 


Process 


BSM 


4300 
4330 
4360 
4370 


gg ^ W+W- ^ tt^ll^2^l2 

gg ^ Zj tlltl 

gg 77 


^ anomalous HVV couplings 



Table 3: Process IDs for gluon induced diboson production at LO (one-loop) QCD. 



2.3 Triple vector boson production in association with a hadronic 
jet 

Finally, the triboson plus jet processes given in Table |4] have also been included at NLO 
QCD leveQfi]. As a process involving multiple electroweak bosons and jets, this is an 
important channel in which to compare experimental data with the predictions of the 
Standard Model. The NLO QCD corrections to the total cross section are sizeable, and 
have a non-trivial phasespace dependence. 



ProcId 


Process 


800 


p^p — )• W~^^^j - 




810 


pp — W 77J - 





Table 4: Process IDs for triboson production in association with a hadronic jet at NLO 
QCD. 



3 New features 

In addition to the new processes described above, several existing procedures have been 
extended and new features added. 

^These processes are disabled by default and must be enabled at compilation using the configure 
option --enable-processes=all or --enable-processes=tribosonjet. 
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3.1 Tensor reduction routines 



The tensor reduction routines for up to 4 external legs have been extended for general 
kinematics. Three and four point tensor integrals are extended to deal with Rank 3 and 
Rank 4 integrals, respectively. Furthermore, one-loop tensor reduction routines for up 
to 6 external legs have been included for the massless case. The tensor integrals are 
implemented following Ref. |5| and can be found in the directory loops/TenRed. 

3.2 Anomalous couplings 

Anomalous gauge boson couplings have been implemented for single vector boson produc- 
tion via VBF, which have been seen to have an effect on some distributions |6|, such as the 
azimuthal angle separation of the tagging jets. The existing diboson process pp — )■ WW 
has been modified to include anomalous VVV couplings, as well as anomalous HVV 
couplings in the gluon-induced contributions. The relevant process IDs are given in Ta- 
ble [Sj To run Vbfnlo with anomalous couplings, the switch ANOM_CPL in the input file 
vbfnlo.dat must be switched to true. The anomalous coupling parameters are then 
input via anomV.dat or (for HVV couplings) anom_HVV.dat. 



ProcId 


Process 


BSM 


120 


{-) 
pp - 


Zjj 


l+l-jj 






121 
130 


(-) 
pp - 

{-) 
pp - 


Zjj 
^W+jj 


vo^ijj 
^^I'ejj 




anomalous couplings 


140 


(-) 
pp - 


-^W-jj 


^'i^ejj 






150 


{-) 
pp - 


-> jjj 








300 


(-) 
pp - 


w+w- 




anomalous VVV and HVV couplings 



Table 5: Process IDs for existing processes which have been extended to include anomalous 
couplings. 



3.3 Spin-2 model 

A spin-2 model has been implemented, using an effective Lagrangian to describe the inter- 
actions of spin-2 particles with electroweak gauge bosons for two cases: an isospin singlet 
spin-2 state and a spin-2 triplet, as described in Ref. |7j. This spin-2 model is imple- 
mented for diboson plus two jets production via vector boson fusion {pp — )■ W~^W~jj, 
pp — )■ ZZjj and pp — )■ W Zjj, process IDs 200 - 230). For these processes a spin-2 reso- 
nance is included in addition to the Standard Model diagrams (i.e. both Higgs boson and 
spin-2 diagrams are calculated). 

The file spin2coupl . dat is used to set the parameters for the spin-2 models. It is 
read if the switch SPIN2 in vbfnlo.dat is set to true, and will only run if the spin-2 
models were enabled at compilation using the configure option --enable-spin2. 
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For the singlet spin-2 field, T'^^, the effective Lagrangian is 

(1) 

and for the spin-2 triplet field, T^" , the effective Lagrangian is given by 

Aripiet = \t,^j {h{D^^)^a^{D^^) + hWi^^B-^) , (2) 

where W and B are the usual electroweak field strength tensors, W and B the dual field 
strength tensors, $ is the Higgs field and is the covariant derivative, fi are variable 
coupling parameters and A is the characteristic energy scale of the new physics. 

In order to preserve unitarity, a formfactor is introduced to multiply the amplitudes. 
The formfactor has the form: 

Here Psp2 is the invariant mass of a virtual s-channel spin-2 particle and 2 are the 
invariant masses of the electroweak bosons. The energy scale Aff and the exponent Uff 
describe the scale of the cutoff and the suppression power. 
The input parameters used by Vbfnlo are 

• F1,F2,F3,F4,F5: Coupling parameters for the spin-2 singlet field. Default values 
are F1=F2=F5=1, F3=F4=0. 

• F6,F7: Coupling parameters for the spin-2 triplet field. Default values are set to 1. 

• LAMBDA: Energy scale of the couplings in GeV. Default value is 1500 GeV. 

• LAMBDAFF : Energy scale of the formfactor in GeV. Default value is 3000 GeV. 

• NFF : Exponent of the formfactor. Default value is 4. 

Note that a graviton corresponds to F1=F2=F5=1 and F3=F4=F6=F7=0. 

Vbfnlo also needs the masses and branching ratios (into SM gauge bosons) of the 
spin-2 particles. 

• SP2MASS : Mass of the spin-2 singlet particle in GeV. Default value is 1000 GeV. 

• MSP2TRIPPM: Mass of charged spin-2 triplet particles in GeV. Default value is 
1000 GeV. 

• MSP2TRIPN : Mass of neutral spin-2 triplet particle in GeV. Default value is 1000 GeV. 

• BRRAT: Branching ratio into SM gauge bosons for spin-2 singlet particle. Default 
value is 1. 

• BRRATTRIPPM: Branching ratio into SM gauge bosons for charged spin-2 triplet 
particles. Default value is 1. 

• BRRATTRIPN : Branching ratio into SM gauge bosons for neutral spin-2 triplet par- 
ticle. Default value is 1. 
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A new process pp — )■ •yyjj (Table [6]) has been added which includes only the resonant 
spin-2 diagrams - this can be compared to the existing process of Higgs boson production 
via VBF, where the Higgs decays into photons (process ID 101). The default values given 
above are intended for processes 200-230. Corresponding values for light spin-2 resonances 
in process 191 can be found in Ref. |7P| 



ProcId 


Process 


BSM 


191 


PP ^ S2JJ ^ 77 


spin-2 resonant production only 



Table 6: Process ID for production of a spin-2 particle S2 with 2 jets via vector boson 
fusion at NLO QCD accuracy. 



3.4 Histograms 

The error calculation for the real emission output can now be controlled via the input 
file histograms.dat. Vbfnlo can calculate the Monte Carlo error for each bin and 
output this to the raw histogram data output for ID and 2D histograms. For the gnuplot 
histogram output only the ID histograms can display the error bar^ 

• CALC_ERROR_GNUPLOT : Enable or disable y-error bars in ID gnuplot histograms. 
Default is false. 

• CALC_ERR0R_1D : Enable or disable y-error bars in raw ID histogram output. Default 
is true. 

• CALC_ERR0R_2D : Enable or disable z-error bars in raw 2D histogram output. Default 
is false. 

Furthermore, Vbfnlo uses a smearing between adjacent bins to avoid artefacts at 
NLO when the real emission kinematics and the corresponding subtraction term fall into 
different bins. As this can lead to remnants at the sharp edges caused by cuts, the 
smearing can be switched off. 

• SMEARING: Enable or disable smearing. Default is true. 

• SMEAR_VALUE: Set the bin fraction where the bin smearing is active. The part that 
is put to the next bin becomes larger when the x-value is closer to a bin border. 
Default is 0.2. 

^Note that in Ref. [t] this process is referred to as process 240. 

'^Error calculation is not implemented for the other (RoOT or TopDrawer) histogram formats. 
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3.5 SUSY options for electroweak corrections 

When running Vbfnlo in the MSSM, it is now possible to set the mass of the Higgs 
bosons in the electroweak loops to either their tree-level value or their corrected value, 
using the flag MH_LOOPS in susy.dat. 

In some areas of the MSSM parameter space, the electroweak loop corrections can 
be the dominant contribution to the cross section. In this case, the squared electroweak 
corrections from the (s)fermion corrections are important and can be included in Vbfnlo 
using the susy . dat flag LOOPSQR_SWITCH. If set to true the amplitude is given by 

s)fermion loop I • (4) 

Note that the loop squared component is only added if |A1(s)fermion loopl is greater than 

10% of |MBorn|. 

4 Other changes 

Since the previous release, VERSION 2.5.0, some changes have been made that alter 
previous results (events, cross sections and distributions). 

4.1 Allowed virtuality of resonance 

In the phasespace generators, the allowed range of the virtuality of a resonance of in- 
termediate vector bosons has been increased. This mainly affects processes where an 
intermediate Z boson decays into a pair of neutrinos - i.e. 

• pp — > Hjj — > ZZjj — > t^l'vVjj via vector boson fusion (process ID 107) and 
gluon fusion (process ID 4107) 

• pp ^ Hjjj — )■ ZZjjj £'^£~i'T'jjj (process ID 117) 

• pp ^ Hjjj ZZjjj (process ID 2107) 

• pp — > ZZjj — > t^l'vVjj (process ID 211) 

This not only affects the cross sections for these processes, but also means that the events 
produced by Vbfnlo-2.6.0 will differ from those produced by Vbfnlo-2.5, even if the 
same random numbers are used. 

4.2 Matrix element H ZZ -^U 

A bug was found and flxed in the implementation of the matrix element calculating the 
decay H ^ ZZ ^ AL 

4.3 Anomalous couplings 

Several changes have been made to the implementation of the anomalous couplings. For 
Higgs production via vector boson fusion (process IDs 100-107) the variable TREEFAC, 
which multiplies the Standard Model contribution to the tree-level HVV couplings, 
has been corrected and altered - now, separate factors for HZZ and HWW are input 
(TREEFACZ and TREEFACW respectively). 
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When working with anomalous HVV couphngs two types of formfactors can be apphed 
which model effective, momentum dependent HVV vertices, motivated by new physics 
entering with a large scale A at loop level. Corrections to the HVV formfactor F2, where 

F2 = -2A'Co{qlqliq, + q2f,A'). (5) 

(where the are the momenta of the vector bosons and Co is the scalar one-loop three 
point function in the notation of Ref. [8|) have been made. The implementation of the 
parametrization described by PARAMETR3 - where the input determining the anomalous 
couplings is in terms of the dimension-6 operators {Ow, Ob, Oww and Obb, see |9|[To]) 
- has also been corrected. 

If anomalous triple (and quartic) gauge boson couplings are being studied, a formfac- 
tor, given by 

F=(l + ^Y\ (6) 



A2 

can be applied in order to preserve unitarity, where A is again the scale of new physics. 
The momentum dependence of an applied formfactor (i.e. s) is now universal for each 
phasespace point, with the invariant mass of the bosons as the scale. This ensures the 
proper cancellations for anomalous contributions and affects both the cross sections and 
the distributions significantly. 

The values of the formfactor scale A and suppression p can be set to different values 
for each input describing the triboson couplings. In the parameterization TRIANOM = 2 
of the L3 Collaboration the formfactor scales for Ak^ and Akz are now separately 
set, and the consistency of related parameters (i.e. Agf , Ak,^ and Anz) is enforced when 
formfactors are applied. 

When processes involving resonant Higgs diagrams (e.g. WWW production) are stud- 
ied with anomalous couplings, the Higgs width is now calculated with the appropriate 
anomalous HVV couplings (the anomalous HVV couplings in the production amplitudes 
were taken into account in previous versions of Vbfnlo). Various corrections have also 
been made to the anomalous triboson couplings in diboson plus jet processes (these were 
incorporated into the intermediate release Vbfnlo 2.5.3). 

4.4 VBF Higgs boson production in association with three jets 

A small bug was found and fixed in the calculation of the processes pp — )■ Hjjj, with 
process IDs 110-117. 
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Appendix A: Process List 



The following is a complete list of all processes available in Vbfnlo, including any Beyond 
the Standard Model (BSM) effects that are implemented. Firstly, the processes that are 
accessed via the vbfnlo executable are given. 



Process 



BSM 



100 


(-) 
pp - 


Hjj 


101 


(-) 
pp - 


Hjj 


102 


(-) 
pp - 


Hjj 


103 


(-) 
pp - 


Hjj 


104 


(-) 
pp - 


Hjj 


105 


(-) 
pp - 


Hjj 


106 


(-) 
pp - 


^Hjj 


107 


(-) 
pp - 


Hjj 



11 JJ 
pt"^^" jj 

r^T~ jj 
bbjj 

w+w- jj ^ i+u.^qu^jj 
zzjj^eiqi+qjj 

ZZjj it^i'^h'^hjj 



> anomalous HVV couplings, MSSM 



(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 



Hjjj 

Hjjj 
Hjjj 
Hjjj 
Hjjj 
Hjjj 
Hjjj 
Hjjj 



11J33 
M"^/^" jjj 

T^T- jjj 

bbjjj 

W+W- jjj it ^h^2^f.2 jjj 

zzjjj^eti^e+qjjj 

ZZjjj tt^i^h^tijjj 



(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 

(-) 
pp 



Hljj 

Hljj llljj 

Hljj /J.+ii'jjj 
Hjjj _^ ^+^-jjj 

H-fjj bb-fjj 

Hjjj _^ W^W~-fjj i't^^i^2^t2ljj 

H^jj ^ zz^jj ^ it^i^t^2in 

Hljj ZZ-ijj ii^ii^e^^e^ljj 
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(-) 








120 


PP - 


-> Zjj 


£+£-jj 


> 


121 


(-) 
PP - 


Zjj 


jj 




130 


(-) 
PP - 


^W+jj 




> 


140 


(-) 
PP - 


■> w- jj 






150 


H 
PP - 









Process 



BSM 



anomalous couplings 



(-) 
PP 



S2jj 17 jj 



only spin-2 resonant production 



H 
PP 

(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 

(-) 
PP 
(-) 



ZZjj iiiiiyt^i^e^jj 
W+Zjj^e+ueJt^2jj 

w-zjj^qv,,itqjj 

W+W+ jj ^ l+ue.itutjj 
W'W- jj £f%£2 % jj 



anomalous couplings, Kaluza-Klein Sz spin-2 models 



> Kaluza-Klein models, spin-2 models 



300 


(-) 
PP - 


■> w+w 




310 


(-) 
PP - 


■> w+z 




320 


(-) 
PP - 


^ w-z 




330 


(-) 
PP - 


^ zz^ 


^2 


340 


(-) 
PP - 


■> W+1- 


^i"i'^i7 


350 


(-) 
PP - 


-> W~i - 


^r%7 


360 


(-) 
PP - 


^ Z-f^ 


^r^^7 


370 


(-) 
PP - 


■> 77 





anomalous and HVV couplings 

anomalous VVV couplings 

anomalous HVV couplings 

anomalous VVV couplings 

anomalous HVV couplings 
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Process 



(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 

(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 
PP 

H 
PP 



Z Z Z y ^-j^ ^3 ^3 

—7- £"'"f£77 
VF~77 —7- ^"^'£77 
Z77 —7- £~^+77 
Z77 f£Z^^77 
777 



BSM 



anomalous couplings, 
Kaluza-Klein models 



(-) 
PP 

(-) 
PP 

(-) 
PP 

(-) 
PP 



W'' 7J ^ J^€7J 

w-zj^qi?,,qi+j 

W+Zj ^ £+z^4£2'4i 



p'p — >■ W'^jjj 
pp — >■ vy 77_7 



> anomalous couplings 
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The processes accessed via the executable ggf lo are given below. 



ProcId 


Process 


BSM 


4100 


(-) 
pp - 


H jj 






MSSM, general 2HDM 


4101 


(-) 
pp - 


Hjj 


77 iJ 




> 




4102 


(-) 
pp - 


Hjj 


^^A*" jj 








4103 


(-) 
pp - 


Hjj 


T^T- jj 








4104 


(-) 
pp - 


Hjj 


bbjj 






> MSSM 


4105 


(-) 
pp - 


Hjj 


w+w- 


jj itiyej2 ^^ijj 






4106 


(-) 
PP - 


Hjj 


ZZjj - 








4107 


(-) 
PP - 


Hjj 


^ZZjj- 


> £fi];ui^i>e2jj 


> 




4300 
4330 
4360 
4370 


99 - 
99 - 
99 - 
99 - 


> W+W 

> ZZ^ 

> Z-f^ 

> 77 






> 


> anomalous HVV couplings 
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